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Abstract 

In a previous paper [1] an Euler angle parameterization for SU(N) was given. Here 
we present a generalized Euler angle parameterization for U(N). The formula for 
the calculation of the volume for U(N), CP N as well as other SU(N) and U(N) 
cosets, normalized to this parameterization, will also be given. In addition, the 
mixed and pure state product measures for iV-dimensional density matrices under 
this parameterization will also be derived. 
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1 Introduction 

Having produced an Euler angle parameterization for SU (N) we now turn our 
attention to explicitly writing down the Euler parameterization for the unitary 
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group, U(N) (which was hinted at in the SU(N) work of [1,2]). Recall that 
U(N) is a subgroup of GL(N, C) = GL N (C), the group of all complex N x N 
matrices with non-vanishing determinant and requiring 2N 2 parameters to 
represent. In this manner we can define SU(N) to be a subgroup of U(N), 
requiring N 2 — 1 parameters to represent, by adding the extra condition that 
any element of SU(N) has unit determinant. We can therefore expect that 
not only will the Euler parameterization of U (N) be easy to produce but also 
the group volume, once we exploit some simple group relationships between 
SU(N) and U(N). 

The importance of such a parameterization and its corresponding volume 
equation, beyond that discussed in [3], is that it gives us the ability to cal- 
culate the measures and volumes for general iV-dimensional pure and mixed 
state density matrices, as well as the volumes of the manifolds of operations on 
pure and mixed states which produce entangled and separable states (which 
are directly related to the volume of separable and entangled states) without 
having to resort to extensive numerical computations as in [4,5]. 



2 Euler Parameterization of U(N) 



The idea behind the parameterization of U(N) is straightforward. Referring 
to our previous work [1], for notations and details, as well as to Nakahara 
[6], Sattinger [7], and others we know the following relationship holds between 
SU(N), U(N) and CP^: 



, SU(N + 1) SU(N + 1) 



CRV = i ■ — = i ! — '—. (2.1) 

U(N) SU{N) x U{1) 1 ' 

The U(l) in the denominator of the second equality is the U(l) element from 
the SU(N + 1) group in the numerator, which we know from [1,2] to be: 

U{1) = U(1) SU(N+1) = e'V+D'-^. (2.2) 

Using the SU(N) parameterization work done previously [1], we can write 
down the Euler parameterization of U (N) quite easily. 

Recall from [1,2,7,9] we know we can write down a semi-direct sum for the Lie 
Algebra for SU (N) as 

L(SU(N)) = L(K)®L(P), (2.3) 
which yields a decomposition of the group, 

V — K • P, (2.4) 
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where V G SU(N). From this work, we also know that L(K) is comprised 
of the generators of the SU(N — 1) subalgebra of SU(N), and therefore K 
will be the U(N — 1) subgroup obtained by exponentiating this subalgebra, 
{Ai, . . . , A(jv-i)2_i}, combined with \n 2 -i and thus can be written as (see 
[1,2,10] for examples) 



K(N - 1) = [SU(N - 1)] • e^ 2 -i* (2.5) 

where [SU(N — 1)] represents the (N — l) 2 — 1 term Euler angle parameteri- 
zation of the SU(N — 1) subgroup. 



We are now ready to look at the U(N) group in general. For a U G U (N) we 
have from equations (2.1) and (2.5) as well as from [1] 

U = K(N) = [SU(N)\ ■ e^V+D 2 -^, (2.6) 

where 



[SU(N)] = ( n A(k,j(N))).( n A(k,j((N-l)))...[A(2,j(2)) 

\2<k<N J \2<k<N-\ J \ i 

X e iA 3 a Af2-(JV-l) . . . e iA (Ar-l)2-l Q Jv2_2 e iA iv2-l a ]v2-l 



= n f n Mkj(m))) 

N>m>2 \2<k<m / 



X 



e «A 3 a iv 2_ (JV _ 1) . . . e * A (iV-l)2-l Q! Jv2_2 e * A ]v2_i a ]v2_i ) 



A(k,j(m)) = e^^-ai+jWe^n^+i 02 !^ 1 



j(m) = < 



m — N, 

E 2(m + 1) m^N. 

0<KN-m-l 



(2.7) 



3 Volume of U(N) 

From [1,11] the volume of SU(N) is known to be 

V SU(N) = 2— tt ^ v/jV J] l^-j. (3.1) 

If we use equation (2.2) then from [1] we can define the following volume for 
U(l) 

SU(N+1)- 
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T T _ T, (A J . -,\ rV (JV + 1)((^ + 1)-1) 

^(w + d = = + !) * y ^(at + i )2+ i 

l2(N + l) . . 

= *V ■ (3 ' 2) 

From equation (2.6) we can write 

V[/(jV) = ^5C/(AT) x ^W)a (jv+1)2 x > (^-3) 
and thus using equations (3.1) and (3.2) we have 



V u(N) =2—n 5 V^I] Tih 1 



fc=i 

AT-l 



fc! / V N 



= 2*^V3V + in(^) (3-4) 

for N > 2. Note that when TV = 1 we generate the volume for the f(l)si/(2) 
group element. Since there can be many different £7(1) 's with different volumes, 
the fact that U(N), when N — 1 gives the SU(2) group element volume 
demands that we limit the use of equation (3.4) to iV > 2. 

With this information in hand, we can now look at the differential volume 
elements, and corresponding volumes of the full range of SU(N) and U(N) 
cosets that are of interest in physics, beginning with the fundamental manifolds 
which define pure and mixed states. 



4 Differential Volume Elements for Pure and Mixed States 



Now that we have an Euler angle parameterization for both SU(N) and U (N), 
for N > 2, we are now in a position to look at the group representations of 
pure and mixed states in terms of our parameterizations. 

In general, the manifold of pure states is given by the sequence of maps: 

U(N - 1) h-> SU(N) i-> CP" -1 . (4.1) 

These are related to the "Grassmannian" manifolds, which are defined as 

cp»- sGW i)= , mN ) .= S , U{N \- (4-2) 

V ' ; 17(1) x U(N — 1) U(N-l) V ; 

On the other hand, the manifold for mixed states (here for rank N density 
matrices with non-degenerate and non-singular eigenvalues) is given by [12,13] 
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M ms = Qn-i x J§^k (4-3) 

where flw-i can be seen as the (N — l)-dimensional solid angle (with ap- 
propriate ranges) derived from the eigenvalues of a suitably parameterized 
(N — 1) -dimensional sphere (see [3,13]), and the factor ([/(l))^ -1 is the max- 
imal torus spanned by the exponentiation of the Cartan subalgebra of the 
group 

([/(l))^ 1 = U(1) SU (2) X U(l) SU{ 3) x • • • X U(1) SU (N) 

= U(1) X3 xU(1) X8 x...xU(1) Xn2i . (4.4) 
One may also notice that M ms is stratified by noting that 
SU(N) 



(U(l)) 



N-l 



CP^" 1 x CP N ' 2 x • • • x CP 1 , (4.5) 



where the x, denotes the (possibly) non-trivial topological product of the 
spaces. These cosets are called flag manifolds and the given topological product 
follows from the fact that the SU (N) groups are products of odd-dimensional 
spheres (see [3] and references within). 

Now, in order to do any "physically" meaningful calculation on either manifold 
we require their measures; measures that can be derived by using the Euler 
angle parameterizations of SU(N) and U(N). It is to this question that we 
now turn our attention to. 



4-1 Pure State Measure 



We know that pure states are in CP^ and from the previous sections that 



)V SU(N + 1) SU(N + 1) 



U(N) SU(N)xU(l) su(N+1) - y - J 

Using the differential volume element for SU (N) from [1] we can immediately 
write down the pure state measure as 

_ dV su{N+1) 



dVsu(N) x dV u[1)su(N+l) 



Ksu(N+i)dat( N+ i)2_i . . . dati 



Ksu{N)da,N 2 -i ■ ■ ■ dai x rfo;( A r +1 )2_ 1 

( J] Ker(k,j(N + l)))da 2N ...da 1 . (4.7) 

\2<fc<Af+l / 
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where from [1] 

f sin(2a 2 ) k = 2, 

Ker(k,j(N + 1)) = cos^-!)) 2 *" 3 an^k-i)) 2<k<N + l, (4.8) 

[cos(o; 2 Ar) sin(o;2Ar) 2Ar_1 k = N + 1, 

with the following ranges 

7T 

< «i < 7r, and < a 2 < — 

IX 

< «2j < ^' - a 2i-l < 27T 

for 2<j<N. (4.9) 

Note that these ranges are from the covering ranges for SU (N + 1) and not 
from SU(N + 1)/Zjv+i which are used to calculate the invariant volume for 
SU(N + 1) (see the appendices in [1] for more details). On the other hand, 
one could use the SU(N + 1)/Z N+1 ranges 

7T 

< a 2j < -, < a 2j -i < 7r, 

forl<j<iV (4.10) 

but then one would need to add a normalization factor of 2 Af_1 in front of the 
product in equation (4.7) in order to generate the correct volume for CP N . 



4-1.1 Example Calculation: Two Qubit Pure State Measure 



It is interesting to note that equation (4.7) for iV = 3 is equivalent to the 
"natural" measure (referred to in [2] ) derived from the Hurwitz parameteriza- 
tion (see [14] and references within). To begin, we define a general vector of a 
random 4-dimensional unitary matrix U (4) as 



Mr))) 



cos(6> 3 ) 
sin(# 3 ) cos(# 2 )e^ 3 



\ 



sm(^ 3 j sin(tf 2 J cos^ije 
\^sin(0 3 ) sin(# 2 ) sin(^ 1 )e^ 1 j 



(4.11) 



where < Q { < tt/2 and < < 2n (i = 1, 2, 3), and rj = {9 h fc}. From this 
vector one can calculate the corresponding Fubini-Study metric (here given as 
in [15]) 

9»v = \(^v + %»), (4.12) 
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where in this case 



'_d_ 



;i4 -!*(»?)> <*(i?)i) 
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■*fa)>, 



(4.13) 



the square root of the determinant of which yields the invariant measure for 
CP 3 under this representation: 



dV ps = Det[ y /g\ 



(4.14) 



= cos(#i) sin(^i) cos (6 2 ) sin (6 2 ) 3 cos(8 3 )sin(83) 5 dO^d^s . . . d9\d4>i. 

The equivalent aspect of our statement comes in when one explicitly evaluates 
equation (4.7) for iV = 3: 



dV r 



dVsu(i) 



ps 



dV SU (3) x dV u(1)su(4) 



= Yi Ker(k, j(4)) \da 6 . . . dct\ 

\2<A;<4 / 

= sin(2o;2) cos^) 3 sin(a 4 ) cos(o; 6 ) sin(a 6 ) 5 <ia; 6 . . . da.\ 

= 2 sin(a 2 ) cos(a 2 ) cos(o;4) 3 sin(o;4) cos(o; 6 ) sin(a: 6 ) 5 <ia: 6 . . . da\ (4.15) 

where the ranges on the ctj's are from equation (4.9). Obviously there's some 
contradiction between this measure and the one given in equation (4.14) but 
any concern it may raise should be eliminated in the following work. 

To begin we note that equation (4.15) can also be derived in the following 
manner that follows the arguments found in [16]. First we define a pure state 
as 

A 



p d = 



A) 





= i(l 4 - v^A 
4 



15 



(4.16) 



\0 I J 

and then apply a U G SU (4) to yield 



P = Up' d U^ = ^(U-V6U\ 15 rf). 
Recalling that a general two qubit density matrix has the form 



(4.17) 



|$(«)> ($(«)| = -(1 4 + V6n-A) 



(4.18) 



we can therefore solve for the components of n and in turn $(ck), via evaluating 
rij = $"l"Aj$ for j = 1, . . . , 15. Doing these calculations yields (dropping an 
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overall multiplicative phase term dependent on the A15 element of the Cartan 
subalgebra found in U(3)) 



|<t>(a)> 



( sin(a 6 ) cos(a 4 ) cos(a 2 )e^( ai+Q3+a5 ) \ 
— sin(o; 6 ) cos(a 4 ) sin(a 2 )e^ ai_a!3 " a5 
— sin(a;6) sin(o;4)e 
cos(a: 6 ) 



(4.19) 



Calculating and taking the determinant of the Fubini-Study metric as before 
but now under this representation yields the following invariant measure for 
CP 3 



dV ps = sin(2a 2 ) sin(o;4) cos(a4) 3 sin(o; 6 ) 5 cos(o; 6 ). 



(4.20) 



One can see that |$(a)) is similar to \^{rj)) but not equal. Therefore in 
comparing the two measures we can only note the following: 

(1) The factor of 2 in equation (4.15) equates to having the range of ol\ run 
from to 2tt rather than its original range set < ol\ < n given in 
equation (4.9) thus allowing one to conceptually equate 0j with a 2 i and 
(pi with a functional form of the a^-iS. 

(2) Equation (4.14) can be generalized to CP^ (see [14]) 

N-l 

dV ps = J] cos(^ fe ) smiekY^dOkdfa (4.21) 
k=i 

which obviously does not have the same form as equation (4.7), but due 
to the invariance of the integral 

TV IT 

r sin(O m cos(£K = r sm(0 cos(O m ^ (4.22) 

JO JO 

does yield the same invariant volume (see the next section and [14]). 

Thus the difference between the two pure state measures is just in the way 
one initially chooses the distribution of the angles 77 and ct in the space CP 3 
(and in CP^ in general). Since we are most concerned with unitary operators 
in SU(N) acting upon pure state density matrices and not within the more 
general U(N) group, we feel that our representation of the pure state measure 
is more useful with regards to the overall Euler parameterization of SU(N) 
and U(N) than the one given in equation (4.21). 
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4-2 Mixed State Product Measure 



From equation (4.3) we can see that, in general, one can write down the 
mixed-state product measure for p = U pdU^ as 

dV ms = dpxd(^j, (4.23) 

where dp defines a measure in the (N — l)-dimensional symplex of eigenvalues 
of p d and d(G/H), where G = SU(N) and H = U(l) su{2) x U(l) su{3) x • • • x 
U(l)su(N)i defines a "truncated" Haar measure which is responsible for the 
choice of eigenvectors of p that ensures dp is "rotational invariant." 

Now as [4,5,14,18,19,20,21] and others have noted, dp is defined via the prob- 
ability distribution induced on the (N — l)-dimensional symplex but there can 
be more than one possible dp that is applicable for a given system since there 
can be more than one usable probability distribution. As Hall noted: 

...[I]f [mixed states] described by density operators are allowed, the require- 
ment of unitary invariance (thus there is no preferred measurement basis 
for extracting information) only implies that the probability measure over 
the set of possible states is a function of the density operator eigenvalue 
spectrum alone. Hence a unique probability measure can be specified only 
via some further principle or restriction, to be motivated on physical or 
conceptual grounds [22]. 



Therefore 



An ensemble of general states of a quantum system is in general described 
by a probability measure over the density operators of the system, given 
that probability measures transform in the same way as volume elements 
under coordinate transformations, and that volume elements are in general 
properties of metric spaces, this suggests that the distribution of density 
operators corresponding to a "minimal knowledge" (i. e. most random en- 
semble of possible states) ensemble may be obtained from the normalized 
volume element induced by some natural metric on the set of density ma- 
trices [22]. 

Thus, the volume measure is defined by the choice of metric, and since the 
metric is invariant under unitary transformations, defining dp comes down 
to determining which metric is the most appropriate in defining a statistical 
distance between two density matrices; especially when one adds the additional 
requirement that the metric satisfy certain criteria for entanglement measures 
(see for example [23,24,25,26,27] and references within). 
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Since there are multiple choices for distance measures between two density 
matrices, and therefore d/j,, and since we want to keep our discussion as gen- 
eral as possible, we shall defer further discussion on d/i to other papers (for 
example those previously cited) and just use the most general form of dfi in the 
spirit equation (4.3) and given in [5,19,20] and references within; the Dirichlet 
distribution: 

dn = rl^ Ai 1 ' 1 ' ' ' A^r'a - £ W-'Ui ■ ■ ■ dA N _, (4.24) 

r(si) • • • r(sjv) p[ 

where X Aj = 1 and 1 > Aj > are just the eigenvalues of p&. The Dirichlet 
distribution provides a means of expressing quantities that vary randomly, 
independent of each other, yet obeying the condition that there sum remains 
fixed. In our case, Sj = s > thus 

dfi = ^4 A i _1 • • • A N 1 dA 1 ...dA N 
JML (s) 

= a s A s r X • • • A^dAi . . . dA N , (4.25) 
where the ranges for the Aj, we conjecture, are equal to 

l>A Ar >-iandO<A 2 ,...,A JV _i<-i. (4.26) 

These ranges disagree with those given by Slater in [28] for = 4 but under 
integration, the difference between those in [28] and ours is a multiplicative 
factor of upon the kernel. The benefit of the ranges given here is that they 
are easily generalized, while those in [28] are not. 

With dfi so defined, we are now free to look at the flag manifold G/H. By using 
equations (2.2) and (2.7) we can see that the coset G/H can be expressed as 

G SU{N) 



H U(l) 

SU{2) 

xU(l) 

5(7(3) X • • • X U(l)su(N) 

U N > m >2 (n 2 <K^(^V)))e ,A%2 - (M ' " • .^V-D'-i^-se^-!^ 

e i\30c N 2_ (N _ 1) . . . e iX (N-l)2-l a N2-2 e iX N 2_ 1 a N 2_ 1 

II ( II MkJ(m))\ (4.27) 

N>m>2 \2<k<m ) 



where A(k,j(m)) is defined in equation (2.7). This coset representation comes 
from the following observation; it allows us to write down the "truncated" 
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Haar measure d(G / H) as 

d 



SU(N) 



U(l) su(2) xU(l) su(3) x...xU(l) 



5(7(3) 



'SU(N), 



SU(N) 



dv um sum x ^(i) st /(3) x • • • x dv u(ihu (N) 

Ksu(N)da N 2_i . . . dati 



where from [1] 



K 



SU(N) 



Ker(k,j(m)) 



o?a7v 2 -(Af-i) • • • da^2_i 
= Ksu(N)da N (N-i) ■ ■ ■ dai 



I] f II Ker(kJ(m))\ 

N>m>2 \2<k<m ) 

'sin(2a 2+j(m) ) 

cos(a 2 (fc-i)+i(m)) 2fe ~ 3 sm(a 2(fc _i) +i ( m) ) 
cos{a 2 ( m -i)+j(m)) sin(a 2 (m-i)+j(m)) 2m ~ 3 



(4.28) 



k = 2, 

2 < k < m, 

k = m, 

(4.29) 



and j(m) is from equation (2.7). 



Now the ranges for the a's can again be either the covering ranges defined for 
the first N(N — 1) a's of the Euler parameterization of SU(N) (see [1]) or the 
SU(N)/Z N ranges 



7T 



< « 2 j < -, < a 2 j-l < 7T, 



for 1 < j < 



N(N-1) 



(4.30) 



which would necessitate adding a normalization factor of 2^ 1 ^ JV 2 ^ 2 to 
Ksu(N) in equation (4.28). 

Depending on which set of ranges are used, a general mixed-state product 
measure can thus be written as 



N-1 



dV m 



: c^r 1 • • • A^r'a - E a^- 1 ^ . . . d\ N ^ 

3=1 

x£ • K S U{N)dati . . . dOLN{N-l), 



(4.31) 



where the Aj are the non-zero eigenvalues of the corresponding A^-dimensional 
diagonal density matrix (see [1,4,5,12,14,28,29] for more details) and £ is the 
necessary normalization constant (equal to 1 if one uses the covering ranges 
for SU(N) and equal to 2 (N - 1 ^ N - 2 ^ 2 if one uses the generic SU(N)/Z N coset 
ranges used in calculating the group volume [1,2]. 
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4-2.1 Example Calculation: Two Qubit Mixed State Product Measure 
For two qubits, equations (4.23), (4.25), and (4.28) yield 

SU{4) \ 



dV ms =dfix d\ 



U(l)sU(2) X U(l)su(3) X U(l)su(4), 



a 



.A'r'Ar'Ar^l - E A*)^ 1 ^ ...^xC- K su{4) da 12 ...da, 



i=i 



a s A s 1 - 1 A s 2 ~ 1 A s 3 - 1 Al- 1 dA 1 ...dA 4 

x £ • sin(2a 2 ) sin(o;4) cos(a4) 3 sin(o; 6 ) 5 cos(o;6) 

x sin(2a 8 ) sin(o; 1 o) 3 cos(«io) sin^a^rfa:^ . . . dot\ 



(4.32) 



where we have used the SU(4) differential volume element from [2] in the 
last step. The ranges of integration for the a,i parameters has already been 
discussed; ideally they should be the covering ranges for S77(4) from [2] so 
that £ = 1. As for the ranges on d/i, recall that for two qubits, pa is given by 
[1,2] 



Pd -- 
where 



/sin 2 (6»i)sin 2 (6»2)sin 2 (e 3 ) ^ 

cos 2 (6»i)sin 2 (6»2)sm 2 (6» 3 ) 

cos 2 (0 2 ) sin 2 (6» 3 ) 



V 















(4.33) 



T<^i<-> cos- 1 (i<^ 2 <-, -<^ 3 <-, 
4 ~ 2' v v^ ; - 2 - 2' 3 - 2' 



thus we have 



l>A 4 >i 0<Ai<- 0<A 2 <i 0<A 3 <i 
4 4 4 4 



(4.34) 



(4.35) 



for the ranges of integration on dp. Notice that one could have also used 
equation (4.26) with AT = 4 to achieve these ranges, but it is instructive to see 
their explicit derivation. 



5 Volume of CP N and SU(N)/(U(1)) 



Now we are in a position to give two different methods for calculating the vol- 
ume for the pure and mixed state manifolds of general A^-dimensional quan- 
tum systems. The pure state manifold volume is quite simple and already 
well known; it is the volume of CP^, while the mixed state manifold volume, 
as we have seen in the previous section, is the product of two different mea- 
sures - one of which is dependent on the initial distribution of states on the 
(AT — l)-dimensional symplex. Therefore in the mixed state case we shall only 
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worry about calculating the volume contribution from the second measure; the 
"truncated" Haar measure, since the volume from the N —1 symplex measure 
can be equated to a general multiplicative constant determined by the initial 
distribution of states on the N — 1 symplex (see for example [28]): 

^Mixcd States = ^Symplex X ^[/(^/([/(l))^- 1 = uV S u (N) / (U (1)) N - 1 ■ (5-1) 

For example, for the two qubit case described previously, a naive calculation 
of uj can be seen to be equal to: 



uj = a 



1 - - - 

/ So So So A r lA r lA r lA 4~ lrfA i rfA 2 rfA 3^A 4 

T{As) / 4- 4 *(-l+4*) \ 
_ 4T(2) V ^ J 

for when s > (note that the integration ranges on A 4 were reversed in order 
to keep uj positive). 



5.1 Volume of CP* 



Using the results for U (N) we can immediately write down the general volume 
for CP^. Using equations (2.1), (3.1) and (3.4) we have 



T/ Vsu(N+l) 

V CP N = 77 

VU(N) 



7T N 



2f^ v ^n n ^(^) 

(5.3) 



We see that this result also comes from the integration of equation (4.7) over 
the ranges given in equation (4.9) or equation (4.10) 



/•••/( II Ker(kJ(N+l))\da 2N ...da 1 = 2 N - 1 7r N ]J V(k,N+l) 

I ranges \2<k<N+l ) 2<k<N+l 

(5.4) 

where V(k, N + 1) is, from [1], 



1 k = 2, 

V( ^ + 1) ^U 2< k <_ N+ , <-» 

Expansion of this product yields 

2<fc<Af+l <± u z,iv z, iv. 
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which when multiplied by the 2 7V ~ 1 7i JV factor gives the volume for CP^ as 
previously calculated. One can also see that integration of equation (4.21) 
using the ranges given for 77 in equation (4.11) will also generate equation 
(5.3) (see [14] for example). 

We should also note the remarkable results: 

00 00 n 

Vol(CP n ) = ]T — = e n « 23.147 

n=0 n=0 n - 

k 

lim J] Vol(CP n ) -> 0. (5.7) 



n=0 



Thus, in terms of our pure state manifold discussion, we can conclude that as 
one increases the dimensionality of the system, there will always be a non-zero 
pure state volume. In the spirit of this result we should also note an interesting 
introduction to the importance of the pure state manifold CP^ for large N, 
especially with regard to quantum entanglement, can be found in [30]. 



5.2 Volume of SU(N)/(U(1)) 



Recall that the action of a group in the adjoint representation produces in- 
teresting orbits; the manifolds of which are called generalized flag manifolds, 
and appear very often in geometric quantization, density matrices, entangled 
states, etc. These manifolds can be represented by the coset SU(N)/U(1) N ~ 1 ; 
obviously then the volume of such manifolds are quite important to our work. 
By using equations (3.2) and (4.4) we can write down the general volume for 
such a coset, SU(N)/U(1) SU{2) x U(l) su{3) x • • • x U(1) SU{N) , as 

y( SU(N) \ _ V SU[N) 



U(l)x 3 x ^(1)a 8 x ••• x U(1) Xn2 J V u(lhs x V U(1K x...V LVK 



- N-l (JV-l)(jV+2) 

2 2 7f 2 



27T * v^TT * ... * Tl^J^Zl 



' L (5. 
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But we know 



N-l 



Thus we can see that 

which is in agreement with [3]. This volume can also be generated via in- 
tegrating the "truncated" Haar measure, using the appropriate ranges and 
normalization conditions, given in equation (4.28) (done in detail in [1]). 

For completeness, and with equation (4.5) in mind, we should note the follow- 
ing: 

7T^-^ = TJ n k . (5.H) 
k=l 

Thus equation (5.10) has the following, equivalent representation 



v ( ?m ) = ^ W (i) = fr 1 + 

{um,, x t/(l) As x ... x U(1), n2 J i\ [aj i\ « 

N-l 

= n ^cp*- ( 5 - 12 ) 

k=l 

So the volume of our flag manifold is nothing more than the product of the 
volumes of the complex projective space CP fc where k < N [3]. Notice also 
that as N increases the volume of the flag manifold approaches, but never 
equals, zero (see equation (5.7)). It is an asymptotic limit which converges to 
zero from the left on M 1 . Thus, since one usually chooses a non-zero probability 
distribution on the N — 1 symplex defining dfi (see equation (4.25)), we can 
conclude that, as in the pure state case, the mixed state volume measure will 
never equal zero unless Vs ymp iex doesl 



6 Other SU(N) and U(N) Coset Volumes 



Beyond the full pure and mixed state manifolds there are numerous other sub- 
manifolds that are of interest in physics; the volumes of which have already 
been calculated (see, for example, [3,11,21,31,32,33] and references within). 
These sub-manifolds and their volumes give us both a way to confirm our 
methodology, as well as offering a systematic, rather than numeric, way of 
computing such quantities. From this work we will then be able to calculate 
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the manifolds that contain the set of entangled and mixed states (either pure 
or mixed) for specific quantum systems [34]. It should be understood though 
that the following volume calculations are specific to the SU(N) and U(N) 
Euler angle parameterization that we have developed and its corresponding 
normalizations via the Cartan subalgebra being used. The general question of 
volume normalization of a manifold, especially when one begins to talk about 
coset manifolds with specific elements of the Cartan subalgebra being removed 
will be the subject of a future paper. 



6.1 Volume of SU(N)/SU(P) x SU(Q) 



To begin, we would like to be able to write down the general volume of the 
SU(N)/SU(P) x SU(Q) coset where N + 1 > P + Q and P, Q ^ 1. To do 
this we can use equation (3.1) to generate 



VSU(N) _ V ', 



W) >< Vsu (Q) 2^ L ^¥^VPXl P k =l{i) x 2^^^^VQUtH h 



(N+1)-(P+Q) N (N+l)-P(P+l)-Q(Q+l) + 2 I N ^-J ( 1 \ C-f 1 , , \-} , , 

= 2 2 * 2 v^SUJS w S* L 

(6.1) 



When N + l= P + Qwe have 



VsU(N) (N+1)-(N+1) (P+Q-l)(P+Q)-P(P+l)-Q(Q + l)+2 

— = 2 2 7T 2 

Vsu(P) X Vs{/(Q) 

Ip + O- 1 P+Q ~ 2 / 1 \ 0-1 
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6.1.1 Example Calculation: Volume ofSU(A)/SU(2) x SU{2) 



Defining N = 4, and P = Q = 2, we get from equation (6.1) the volume of 
the coset SU{A)/SU{2) x SU(2) 



VSU(4) ~$=± 4(5)-2(3)-2(3) + 2 / 4 VrV 1 X ~" ~ 



^ ■ j^n j n«n« 



Vs!/(2) x Vsup) " V( 2 )( 2 )"A* ! 
= \f2 *7t 5 *^|*1!>k1! 

7T 5 



k=l k=l 



QV2' 



(6.3) 



Which is equivalent to the volume of SU(A), \^2ir 9 /3, divided by the square 
of the volume of SU(2), 2n 2 , as expected. It should also be noted that this is 
the volume of the manifold that is comprised of all non-local transformations 
which can be implemented on a two qubit system. 



6.2 Volume ofSU(N)/U(P) x 17(1) 



Beyond the general volume of CP , general flag manifold, and the previous 
SU(N) coset, we would like to be able to write down the general volume of 
the SU(N)/U(P) x U(l) coset where N - 1 > P + 1 and P ^ 1. To do this 
we can use equations (3.1) and (3.4) as follows 



N-l (JV-l)(iV+2) rj-r N _ ± 

VSU(N) 



k\ I 

J—. (6.4) 
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The problem we now face is how to define U(l). If we use equation (3.2), here 
now defined for SU(N), we would generate 



SU(N) 



N-l (N-l)(N+2) r— N _ x 

2 2 7T 2 yN\Yk=\ 



fc! 



(jy-l)-(P + l) (iy2 +J V-2)-(p2 + p +2) IN — I ( \\ P ^ 

2 2 * ' VpttS U S«- 



fc=i v-/ fc=i 

(6.5) 



If we demand that TV — 1 = P + 1, we can simplify the product terms 



v-i / i \ P-l P+l / I \ P-l 

sWn^n( M )n« 

1!*2!*---*(P-2)!*(P- 1)! 



1! * 2! * • • • * (P - 2)! * (P - 1)! * P * (P + 1)! 

' (6-6) 



P!(P+ 1)! 

as well as the powers and other factors. Therefore, for this case we have 



(JV-1)-(P+1) (iV 2 +JV-2)-(P 2 + P+2) 



Vsu(N) _ 2 2 7T 2 N-l 



Vu(P)xV u{1)su(N) P!(P + 1)! VP+1 

r 27V-3 



7T 



(N-2)\(N- 1)! 

7T 2P+1 



P!(P + 1)!' 



(6.7) 



Depending on which parameter is used. 

Now, if in using equation (3.2), we now define U(l) for SU{M), M < N, we 
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would generate 



SU(N) 



Vu(P) x V u(1)smil) 



7T 



\M(P + l)l\{k>)l\ *- 

(6.8) 



and if we demand N — 1 — P + 1, we can simplify, yielding 



(JV-1)-(P+1)) (iV 2 + iV-2)-(P 2 + P+2) 
2 2 7T 2 



P!(P + 1)! 

2N-3 



AT(M - 1) 
^ M(P+1) 



7T 



iV(M-l) 
(JV-2)!(JV- 1)!^ M(iV- 1)' 



(6.9) 



which reduces to equation (6.7) when M = N. Therefore, depending on which 
U(l) we use, we will generate a different volume; the ratio between any two 
being equal to 



Vsu(n)/Vu(p) x Vu(i )suix) Y(X - 1) 
Vsu { n)/V u{p) xV u{1)su(y) \X(Y-1Y 



(6.10) 



6.2.1 Example Calculation: Volumes of SU(4)/U(2) x U(l) S u(i) fori = 2,3,4 

Defining N = 4 and P = 2 (thus satisfying N — l = P + l we get from equation 
(6.9) the volume of the coset SU(A)/U(2) x U(l) SU (i) when % = 2 



V, 



5(7(4) 



7T 



2*4-3 



V[7(2) X V U(1)su(2) 

when i = 3 

Kst/(4) 



4(2-1) 7T 5 /2 



7T 



(4-2)!(4- 1)!\ 2(4- 1) 12 V 3 6^6' 



(6.11) 



7T 



2*4-3 



Vu(2) X Vu(l) 



SU(3) 



4(3-1) 7T 5 



7T 



(4-2)!(4- 1)!\ 3(4- 1) 12 V 9 9^2' 



(6.12) 



and when % — 4 the volume of the coset SU{A)/U{2) x U(l)su(4), using equation 
(6.7) now, is 



VsU(4) 



2*4-3 



7T 



Vu(2) x V C7( i )s[7(4) 



(4-2)1(4-1)! 12' 



(6.13) 
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6.3 Volume ofSU(N)/U(P) x U(Q) 



Now we would like to be able to write down the general volume of the SU (N) / U (P) x 
U(Q) coset for iV — 1 > P + Q and P, Q ^ 1. To do this we can use equations 
(3.1) and (3.4) as follows 



V, 



SU(N) 



Vu iP) x V u(m 2 f ^^p+rn^i Ji j x 2%^VQTTYltl(i) 

(6.14) 

Simplification yields 

2 2 7T 2 [[ k=1 



v. 



SU(N) 



k\ 



Vu(p) X Vlf(Q) P+Q P(P+1)+Q(Q + 1) 

V ' 2 2 7f 2 



v / (p + i)(Q + i)nL- 1 1 nK (±) 



(JV-1)-(P+Q) (Af-l)(iV+2)-P(P+l)-Q(Q + l) / iV 
= 2 2 7]" 2 



(P + 1)(Q + 1) 



W-l / \ P-l Q-l 



x n ky n^n^ (6.15) 

k=l VW fc=l k=l 

For the special case when TV — 1 = P+Q we can go further and eliminate the N 
dependence in the above volume, thus yielding (in one possible representation) 



Vsu(N) _ ^(p+q+pq) / P + Q + l tj- ( 1 



P+Q / i \ P-l Q-l 



Vu(p) x V V (Q) 



™ri)n( n jn«n« 



_ ^(p+q+pq) / -P + Q + i tt 1 1 



fc=l k=l k=l 

P+Q / i \ Q-l 



O.i Example Calculation: Volume of SU(9)/U(A) x £7(4) 

Defining N = 9, and P = Q = 4, thus satisfying iV — 1 = P + Q, we get from 
equation (6.16) the volume of the coset SU(9)/U(A) x £7(4) to be equal to 



4+4 / i \ 4-1 



VW(9) r 4+ 4 + 4*4) 4 + 4 + 1 t^/l 



7T 



V(4 + l)(4 + l) fc 1 = 4V^//. 

7T 24 

= 58525286400000' 
which is what one would get if they used equations (3.1) and (3.4) separately. 
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6.4 Volume ofSU(N)/UU U{Pi) x Yl v j= i U(l) su{Zj 



We are now ready to write down the volume for the most general of cosets 
that we are interested in, SU(N)/H* =1 U(Pi) x H%iU(i) su ^ z ^ where 

£^ + £i = £^ + y< N-i, Pi^i, (6.18) 

i=l j=l i=l 

and 

^(1 W,j e {^(l) W ( 2 ), ^(1 W ■ ■ ■ > ^(1)^}. ( 6 - 19 ) 

where there is no necessary order in the sequential choice of U(l) su ^ z y 
To begin we note the following using equation (3.4) 

V(f[U(P z )) = f[V u{Pt) 



1=1 



fi at: — 7 ^i-r 1 ( 1 * 



^.^n(^TTn(3). (-) 

We also can simplify the second volume of the three we need via generalizing 
equation (3.2) 



y y 



2Z j 



v(U u(i)s U{Zj) ) = n vu(i )mnzj) = n -v z , 

j=l j=l j=l V 3 



y 



We are now in a position to write down the volume for SU(N)/ nf=i f(-Pi) x 
nj=i ^(l)st/(z-)- Using equations (3.1), (6.18), (6.20), and (6.21) we have 



/ SU(N) \ N -( 1+ V+T,i=l P ^ (N-l)(N+2)-(2y+J2 X i ^ 1 + 

^ „„ =2 2 yf 2 



.nti^xn^i^iw, 



xVN -. V ' T^\ - ( 6 - 22 ) 



uti [VK+Tim 



Pi-i i j_ 

fc! 



For the special case when the "<" in equation (6.18) is replaced by "=" we 
have 
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/ SU(N) \ (N-l)(N+2)-(2y+Y^ X i _ 1 Pi(P i +l)) 

^ — (6 - 23) 



n^U 1 ( l\ j nj=i 

xv 7 ^ 7 V 7 



n^i v^+TnjgTiMir 



One could continue simplifying equation (6.22) but it would be only worth- 
while if additional knowledge concerning Zj and Pi was available. 



6.5 Grassmann Volume 



The general Grassmann manifolds, of which CP^ is a special case (see equation 
(4.2)), have the following definition for N > M 



G{N ^= U(MuZ-MY (6 - 24) 



Using equation (3.4) we can write down the general expression for the volume 
of almost any Grassmann manifold: 



T/ V U(N) 
VG(N,M) ~~ 



Vjj(M) X Vu(N-M) 



2«^VivTTnt 1 1 (V 



N-l / -i \ M-l jV-M-1 



= 7T 



M(N-M) / iV" + 1 

V (M + 1)(7V-M + 



W-l / 1 \ M-l N-M-l 

i)n(i)n« n « 



= 7T 



"klra n m n «■ c*) 



k=M \ / k=l 
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The reason for the "almost" above is that for M = 1 we do not regain the 
volume for CP N_1 that we originally calculated in equation (5.3): 



v CPN -. = v am) = ^y (1+1)(iV _ 1+1) n 



N + l J ^V1_ 

jfe! 



N-l / i \ iV-1-1 



JV- 



— n 

rN-1 



l N + l"' 1 1 



k=i 

N-l / -i \ N-2 




- n e n« 



iV+ 1 




(iV-1)!- 



(6.26) 



We are "off" by a factor of y (-/V + 1)/2N which occurs because of the following 
reason: equation (3.4), for N — 1, yields 27r which is correct if one is looking 
for the volume of the SU(2) variant of U(l) (see equation (3.2) for N = 1), 
but that is not the case for the U(l) components of greater SU(N) groups 
(again see equation (3.2) for N > 2), which is the case here. In equation (6.26) 
we get the factor of 2tt from the U(l) component, but without the additional 
contraction term due to the \n 2 -i Cartan subalgebra component of U(N) 
from which the U(l) term is defined! 

The flaw in equation (6.26) can also be seen from equation (4.2): 
cp n-i = G(N, 1 



v-j _ r<, s 1 , U ( N ) SU(N) x U(l) SU(N) 



17(1) x U(N - 1) 17(1) x U(N - 1) U(N-l)' 

(6.27) 

The U(l) term in the numerator is the same as the Z7 (1) term in the denomina- 
tor and as such, in any representation, cancels out, thus leaving the standard 
coset relationship for CP W_1 which, from equation (5.3), does yield the correct 
volume for CP^ -1 . Therefore, if we use equation (3.2) for the U(l) term in 
equation (6.25) (combined with equation (3.4) for the other two terms) when 
M — 1, we will generate the correct volume for CP^ 1 



Vg(n,i) 



VjJ{N) 



2f^ v ^nnf7 1 1 (^ 



V U{1) x V U{N - 1} ^y^ry x 2 E f 1 n iE ^ m VNn^ 1 2 ( £ 



(N-1)V 



(6.28) 



Therefore, in general, if we demand that M ^ 1 then equation (6.25) will 
correctly produce the Grassmann manifold volumes (see [3] and references 
within). 
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7 Conclusion 



Using the volume equations given herein, we are now in a position to explic- 
itly write down the measures and volumes for the whole range of manifolds 
which occur in discussions concerning separability and entanglement of multi- 
particle systems. Therefore, this work allows us to explicitly write down the 
volumes of the manifolds of the local orbits of a given state \ip) with respect to 
some transformation U G SU(N) (or more generally U(N)), in a manner that 
we hypothesize also elucidates the topology of the manifolds as well [3,33]. 
Applications beyond quantum information theory are also possible [3]. 
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